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ABSTRACT 

We present a Monte-Carlo event generator for simulating chargino pair- 
production at the International Linear Collider (ILC) at next-to-leading order 
in the electroweak couplings. By properly resumming photons in the soft and 
collinear regions, we avoid negative event weights, so the program can sim- 
ulate physical (unweighted) event samples. Photons are explicitly generated 
throughout the range where they can be experimentally resolved. Inspecting 
the dependence on the cutoffs separating the soft and collinear regions, we 
evaluate the systematic errors due to soft and collinear approximations. In 
the resummation approach, the residual uncertainty can be brought down to 
the per-mil level, coinciding with the expected statistical uncertainty at the 
ILC. 
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1 Introduction 



The MSSM, the minimal super symmetric (SUSY) extension of the Standard Model (SM), is 
a promising candidate for a theory of electroweak interactions pQ. In this model, the Higgs 
sector is stabilized against power divergences in radiative corrections, proton stability suggests 
a discrete symmetry that provides us with a dark-matter particle, and the renormalization- 
group evolution of couplings is precisely consistent with gauge-coupling unification (GUT) at 
an energy scale of the order 10 16 GeV. 

A solid prediction of the MSSM is the existence of charginos X11X21 the superpartners of 
the W ± and the charged-Higgs H bosons. In GUT models their masses tend to be near the 
lower edge of the superpartner spectrum, since the absence of strong interactions precludes large 
positive renormalizations of their effective masses. Thus, if any superpartners are accessible 
in e + e~ collisions at a first-phase ILC with cm. energy of 500 GeV, the lighter chargino xf 
is likely to be pair-produced with a sizable cross section. In many models, including popular 
supergravity-inspired scenarios such as SPSla/SPSla' the second chargino xf a ^ so be 
accessible at the ILC, at least if the cm. energy is increased to about 1 TeV. Similar arguments 
hold for the neutralinos, the superpartners of neutral gauge and Higgs bosons. The lightest 
neutralino is possibly the lightest superpartner (LSP) and therefore the dark-matter particle 
present in the MSSM. 

The precise measurement of the chargino parameters (masses, mixing of Xi with xf > an d 
couplings) is a key for uncovering any of the fundamental properties of the MSSM that we 
have mentioned above. These values give a handle for proving supersymmetry in the Higgs 
and gauge-boson sector and thus the cancellation of power divergences. Charginos decay either 
directly or via short cascades into the LSP, and a precise knowledge of masses and mixing 
parameters in the chargino/neutralino sector is the most important ingredient for predicting the 
dark-matter content of the universe. Finally, the high-scale evolution of their mass parameters 
should point to a particular supersymmetry-breaking scenario, if the context of a GUT model 
is assumed (cf. |3J). In all these knowledge of parameters with at least percent-level 

accuracy is necessary. 

At the LHC, this is a difficult task since charginos are accessible mainly in complicated decay 
cascades of colored superpartners with substantial background, while direct pair-production is 
suppressed jl]. The ILC provides much cleaner production channels and decay signatures with 
low background, so the required precision will be available at the ILC |5|b'j . To match this 
experimental accuracy, there is obvious need for theoretical predictions with next-to-leading 
order (NLO) accuracy in the electroweak couplings. The predictions have to be implemented 
in the simulation tools that are actually used in the experimental analyses (e.g. see [Jj). 

At leading order (LO), chargino production at the ILC is given by the tree-level diagrams in 
Fig. Q and events can be generated using the narrow- width approximation where all processes 
are factorized in on-shell 2 —>■ 2 production and a cascade of on-shell 1 — ► n decay processes. 
The helicity amplitudes can be expressed in analytic form (cf. jH]), and the process is available 
in various computer codes 

The NLO corrections include 1 (i) loop corrections to the SUSY production and decay pro- 

1 This describes the multiple-pole approximation JOj; recent complete NLO calculations in the SM have 
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Figure 1: Feynman graphs for chargino pair production at the ILC. 





Mass 


Width 


xt 


183.7 GeV 


0.077 GeV 


xi 


415.4 GeV 


3.1 GeV 



Table 1: Chargino masses and widths for the SUSY parameter set SPSla'. 

cesses, (ii) nonfactorizable, but maximally resonant photon exchange between production and 
decay, (hi) real radiation of photons, (iv) off-shell kinematics for the signal process, (v) ir- 
reducible background from all other multi-particle SUSY processes, and (vi) reducible, but 
experimentally indistinguishable background from Standard-Model (SM) processes. So far, no 
calculation and simulation code provides all NLO pieces for a process involving SUSY particles. 

In Ref . [E] , three computer codes have been presented and verified against each other that 
simulate off-shell multi-particle processes at tree-level, both for the SM and the MSSM. As 
generators of unweighted SUSY event samples, they thus cover (iv), (v), and (vi). In particular, 
the program described in this paper is implemented as an extension to the WHIZARD event 
generator With beamstrahlung, resummed initial-state radiation, arbitrary polarization 
modes and standard parton-shower and hadronization interfaces being included, this generator 
is well suited for ILC physics studies. 

In this paper, we describe the extension of the tree-level simulation of chargino production 
at the ILC by radiative corrections to the on-shell process, i.e., we consider (i) in the above 
list and consistently include real photon radiation (iii). This is actually a useful approximation 
since in many MSSM scenarios charginos, in particular Xi , are quite narrow (cf. Table QJ, so 
nonfactorizable contributions are significantly suppressed and decay corrections can be sepa- 
rated from the corrections to the production process. 

We emphasize that for the simulation of physical (i.e., unweighted) event samples, it is 
essential that the effective matrix elements are positive semidefinite over the whole accessible 
phase space. The QED part of radiative corrections does not meet this requirement in some 
phase space regions. Methods for dealing with this problem have been developed in the LEP1 
era [T4"| . While these methods are also applicable for the ILC situation, they need a thorough 
reconsideration since the ILC precision actually exceeds the one achieved in LEP experiments. 

explicitly verified the validity of this approximation in the signal region. 
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Figure 2: Chargino pair production at the ILC: Dependence of the differential distribution in 
polar angle cos 9 between e~ and X\ f or different helicity combinations. The labels indicate Xi 
and xf helicity; the electron/positron helicity is fixed to — K 

2 Fixed-Order Simulation of Chargino Production 
2.1 Lowest Order 

In the MSSM, the charginos xf,xf are mixtures of weak gauginos w ± and higgsinos h ± . 
The production processes in e + e~ collisions are thus connected by SUSY transformations to 
e + e~ — > W + W~ and e + e~ — > H + H~; the contributing Feynman diagrams consist of s-channel 
e + e~ annihilation via Z and photon and t-channel exchange of an electron-sneutrino (Fig. 
Since no massless particles are exchanged in the t-channel, the electron mass can be neglected at 
tree level throughout the phase space. The square of the absolute value of the matrix element, 
integrated over the phase space T which is parameterized by production angles 9, 0, defines the 
Born cross section <TB orn : 



We suppress the dependence on particle masses M^, M z , etc. 

At the ILC, the possibility of polarizing electrons and positrons in the initial state will allow 
for separately measuring individual (squared) helicity amplitudes. Selecting a standard MSSM 
parameter point SPSla' (cf. App. O and Table and a collider energy of 1 TeV, in Fig. 121 
we display the angular dependence of the cross section for the dominant helicity combinations 
with v exchange in the t channel. As the amplitudes are oc (1 ± cos 8), sin 9 respectively, they 
can become zero for 9 = ±7r, 0. 




(1) 
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2.2 NLO Corrections 



The one-loop corrections to the process e~e + — ► XiXj with j — 1,2 have been computed in 
Ref. pH], using the FeynArts/FormCalc package ^Hj for the evaluation of one-loop Feynman 
diagrams in the MSSM. An independent calculation with consistent numerical results has been 
presented in [T7] . These calculations include the complete set of virtual diagrams contributing 
to the process with both SM and SUSY particles in the loop. The collinear singularity for 
photon radiation off the incoming electron and positron is regulated by the finite electron 
mass m e . As an infrared regulator, the calculation introduces a fictitious photon mass A. The 
interference of these diagrams with the Born term defines the 'virtual' contribution 

a virt (s, X 2 ,m 2 e )= J dT 2 [2Re (M Born (s)*Mi -loop "J? 

X 2 ,m 2 e ))}. (2) 

The dependence on the fictitious parameter A is eliminated by neglecting contributions pro- 
portional to powers of A and adding real photon radiation with energy E 1 < AEj, where E 1 is 
defined in some reference frame, usually the cm. frame. Hence, the residual logarithmic depen- 
dence on A is cancelled in favor of a logarithmic dependence on AE 1 . This correction can be ex- 
pressed as a universal factor / so ft(— r 1 ) (EH)- The 'virtual + soft' contribution a v+s (s, AE 7 , ml) 
is thus given by 

/soft(^) \M Bo m(s)\ 2 +2Re(M Born l-loop 

(s,X 2 ,rnl))]. (3) 

In a real experiment, there is always a finite energy resolution AE^ xp for photons, and combining 
soft and virtual photons below this cutoff is justified. For the simulation one would choose 
AE 1 < AE^p. 

This result is complemented by the 'hard' contribution a" 2 -> 3( s > AE 7 , m 2 ,), i.e., the real- 
radiation process e~e + — > XiXjl integrated over photon phase space down to the energy 
resolution AEy. 

a 2 ^ 3 (s,AE^m 2 e ) = [ dT 3 \M 2 -* 3 {s, m\)\ 2 . (4) 

Jae 1 

The sum, which can be expressed as a total cross section or, e.g., as a differential distribution 
in the chargino polar angle 9, should not depend on the photon-energy cutoff: 

a tot (s,m 2 e ) = cr Bo rn(s) + Ov+s(s, AE 1) ml) +a 2 ^ 3 (s,A£ 7 ,m^) (5) 

However, the dependence on AE 1 cancels only approximately since positive powers of AE^ are 
neglected in the v+s term but not in the 2—^3 process. 

2.3 Collinear Photons 

While photons with large energy and large angle can be experimentally resolved and must be 
explicitly generated by the Monte-Carlo simulation program, photons collinear to the incoming 
electrons cannot be detected. (Since the outgoing charginos have substantial mass, a collinear 



a v+s (s, AE^,m 2 e ) = / dT 2 
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approximation for final-state radiation is not needed.) As usual, we break down the (hard) 
2^3 cross section into a collinear and a non-collinear part, separated at a photon acollinearity 
angle A9 1 relative to the incoming electron or positron: 

, TTl e ) C~hard, non-coil 

(s, A£ 7 , A0 7 ) + a hardjCO ii(s, A£ 7 , A6» 7 , m 2 e ), (6) 
where in the non-collinear part the electron mass can be neglected. The last term is approx- 

2 

imated by convoluting the Born cross section with a structure function f(x;A9 J ,^ s -) : with 
x = 1 — 2Ey/ \fs being the energy fraction of the electron after radiation, 

0Wd,coii(s, AS 7 , A0 7 , m 2 ) = / dT 3 \M 2 ^ 3 {s,m 2 e )\ 2 

= £° dxf(x;A9 7 ,^) J dT 2 \M B oUxs,m 2 e )\ 2 . (7) 

2 

The structure function f(x; A# 7 , — a ) contains two pieces / + , /~ (|21I23|) that correspond to he- 
licity conservation and helicity flip, respectively; each one is convoluted with the corresponding 
matrix element. The cutoff AE^ is replaced by xq = 1 — 2 AE y / \fs. In this approximation, 
positive powers of A9 1 are neglected. 



2.4 Simulation 

Combining the above, the cross section is given by 

a tot (s,m 2 e ) = J dxf eS (x u x 2 ;AE 7 ,A9 1 ,^) J dT 2 \M cS (s,x 1 , x 2 \m 2 e )\ 2 

+ [ dT 3 \M 2 ^ 3 (s)\ 2 , (8) 

where we define 

/eff(xi, x 2 - AEj, Afl 7 , =£) = 5(1 - zi) 5(1 - x 2 ) 

+ «J(l-x 1 )/(x 2 ;A0 7S ^)0(x o -x 2 ) 

+ /(x 1 ;A0 7 ,^)5(l-x 2 )^(x o -x 1 ) (9) 

and 

|/W cff (s,Xi, x 2 ;m 2 )| 2 = [l + / soft (A£ 7 ,A 2 )#(xi,:r 2 ))] |-M B orn(s)| 2 

+ 2Re [7W B orn(s) -Mi_ioop(s, A 2 , m 2 )] flfo, x 2 ) (10) 

with 9(xi,x 2 ) = 9(xi — xq) 9{x 2 — x ). 

This structure is suitable for implementing it into an event generator. In WHIZARD, for 
instance, there is an interface for arbitrary structure functions f(xi, x 2 ) that can be convoluted 
with the Born squared matrix element. We insert the above effective radiator function f e g as 
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a 'user-defined' structure function and replace the Born matrix element as computed by the 
matrix-element generator, O'Mega ^H], by the effective matrix element defined above. The 
latter is computed by a call to the FormCalc-generated routine. 

In order to account for the ^-function part contained in the radiator function, for the Monte- 
Carlo sampling of x values the Xi range is separated into two regions each, one for Xi < xq and 
the other one for Xi > xq. For each Xi, the first region is mapped such as to maximize the 
efficiency of event generation. If the sampled point ends up in the second region, x^ is set equal 
to 1 before the matrix element is evaluated as demanded by the S function. The relative weight 
of the two regions is given by 

w(x > Xq) : w(x < Xq) = 1 : / dx f(x; A0 7 , ^j-). (11) 

Jo 

For a consistent first-order calculation, we have to avoid the radiation of two (collinear) photons. 
Therefore, the radiator function f e g is zero in the region with x\ < Xq and x 2 < x , and in the 
2^3 process, no convolution with structure functions is applied. 

Implementing this algorithm in WHIZARD, we construct an unweighted event generator. With 
separate runs for the 2 — ► 2 and 2^3 parts, the program first adapts the phase space sampling 
and calculates a precise estimate of the cross section. The built-in routines apply event rejection 
based on the effective weight and thus generate unweighted event samples. 

For the 2 — > 2 part convoluted with a structure function, WHIZARD can optionally represent 
the missing collinear energy by a real photon in the event, with pt generated according to 
the correct logarithmic distribution up to the cutoff angle A9 1 . Thus, if there is any energy 
available for radiation, the actual events contain a photon in addition to the chargino pair 
regardless whether the event has been generated in the 2 — > 2 or 2 — > 3 part. 

On the technical side, for the actual implementation we have carefully checked that all 
physical parameters and, in particular, the definition of helicity states are correctly matched 
between the conventions ^H] used by O'Mega and WHIZARD [20 , and those used by FormCalc 
(cf. e.g. EH)- 

2.5 Where this Approach Fails 

Numerically, the modified WHIZARD code reproduces the total cross section at fixed next-to- 
leading order in a as presented in Ref. [T^]. In principle, this makes the NLO result available 
for physics simulation. However, in the soft-photon region the fixed-order approach runs into 
the well-known problem of negative event weights |22f23j . 

While for any fixed helicity combination and chargino scattering angle the differential cross 
section is positive if we include the virtual contribution and integrate real soft photons up to 
a finite cutoff AE^, the 2 — > 2 part of the NLO-corrected squared matrix element is positive 
definite by itself only if AE^ is sufficiently large, say — j3- = 1CT 2 (i.e., AE 1 = 10 GeV for 

t/s = 1 TeV). If we lower the cutoff, say < 10~ 3 , for each helicity combination the effective 
2^2 matrix element becomes negative within some range of scattering angle, compensating 
the 2^3 squared matrix element that for such a small cutoff overshoots the total NLO 
differential cross section. This is illustrated in Fig. El 
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Figure 3: Effective squared matrix element (arbitrary units) for e~e + — > xlx^ as a function of 
the polar scattering angle 9 at y/s = 1 TeV. Left figure: Helicity combination — h H — ; right 

figure: h ++. Solid line: Born term; dashed: including virtual and soft contributions for 

0.5 GeV. 



AE 7 = 10 GeV; dash- dotted: same with AE^ 



However, the experimental resolution can be better than that, and for large AE^ power 
corrections to the soft approximation are not completely negligible; therefore we are interested 
in letting AE^ — > 0. To obtain unweighted event samples in the presence of negative- weight 
events, an ad-hoc approach is to simply drop such events before proceeding further. The 
numerical consequences of such a truncation are discussed below. 2 

Instead of slicing phase space and introducing a separate treatment of the soft and collinear 
regions, an alternative approach uses subtractions in the integrand to eliminate the singularities 
before integration [24j; the subtracted pieces are integrated analytically and added back or 
cancelled against each other where possible. This method does not require technical cutoffs and 
therefore allows to get rid of cutoff-induced artefacts in the result. Unfortunately, the subtracted 
integrands do not necessarily satisfy positivity conditions either, so the transformation into an 
unweighted Monte-Carlo generator is not straightforward. (Of course, the method may be used 
to construct a weighted event generator.) We do not consider this method in the present paper. 



3 Resumming Photons 
3.1 Leading Logarithms 

The shortcomings of the fixed-order approach described above are associated with the soft- 

E 2 

collinear region E 1 < AE 7 , # 7 < A# 7 , where the appearance of double logarithms ^-ln-^ln# 7 

2 Strictly speaking, events with a weight that can become negative do not preclude detector simulation and 
further analysis, but the event samples containing them are not a possible outcome of a physical experiment. 
In particular, simulating the detector response to such events is a waste of resources. If a method is available 
that eliminates them, it is the preferable choice. 
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invalidates the perturbative series. However, in that region higher-order radiation can be re- 
summed [213 ■ The exponentiated structure function /isr [26] given in Eq. (J2*lj) that resums 
initial-state radiation, 




includes all-order photon radiation in the soft regime at leading-logarithmic approximation 
and, simultaneously, correctly describes collinear radiation of up to three photons in the hard 
regime. It does not account for the helicity-flip part /~ (|2l5|) of the fixed-order structure 
function; this may either be added separately or just be dropped since it is subleading in the 
leading- logarithmic approximat ion . 

In this description of the collinear region, there is no explicit cutoff AE^ involved, and 
collinear virtual photons connected to at least one incoming particle are included. The latter 
part is effectively smeared over small photon energies, such that the +- distribution singularity 
of the finite-order result is replaced by a power-like behavior with a finite limit for x —*■ 1. 
Stated differently, the cancellation of infrared singularities between virtual and real corrections 
is built-in (for collinear photons), so that the main source of negative event weights is eliminated. 



3.2 Matching with NLO 



We combine the ISR-resummed LO result with the additional NLO contributions described 
in the previous section. To achieve this, we first subtract from the effective squared matrix 
element, for each incoming particle, the contribution of one soft photon that is contained in 
the ISR structure function (and has already been accounted for in the soft-photon factor), 



/soft,iSR(A-E 7 , A# 7 , m e 



dx 



■f'o 



1 



X 



7 ( 2 ln(l - x ) + x + ~x\ 



l-x J + 4 V v ' 2 
Here, r\ is defined in Eq. (J22J), and the +- distribution is represented, e.g., by 



g(x) 



lim 



0(1 — x — e) g(x) — 5(1 — x — e) 



1-e 



g(y) dy 



(13) 



(14) 



After this subtraction we have 



\M 



v $\s\ AE 7 , A0 7 ,m; 



2\|2 



l + Zsoft^' 



A ,-2/ softiIS R(A£ 7 ,A0 7 ,^)j \M Bovn (s)\ 2 

+ 2Re [A^Born(s) A*l-loqp(s, ^ , ™e)] , (15) 



with s being the cm. energy after radiation. This expression contains the Born term, the 
virtual and soft collinear contribution with the leading-logarithmic part of virtual photons and 
soft collinear emission removed, and soft non-collinear radiation of one photon; it still depends 
of the cutoff AE T Convoluting this with the resummed ISR structure function, 

^v+ s ,isr(s, A£ 7 , A6> 7 , m 2 e ) 

= [ dxi/isROn;A0 7 ,^) / dx 2 / ISR (x 2 ;A0 7 ,^) / dT 2 \M eS (s; AE 7 , A0 7 , m 2 e )\\ (16) 



S 




Figure 4: Effective squared matrix element (arbitrary units) with the one-photon ISR part 
subtracted, for e"e + — ► X1X1 as a function of the polar scattering angle 6 at y/s = 1 TeV. Left 

figure: Helicity combination h H — / right figure: h ++. Solid line: Born term; dashed: 

including virtual and soft contributions for AE 1 = 10 GeV; dotted: same with AE^ = 0.5 GeV. 
The collinear cutoff is fixed at A# 7 = 1°. 

we obtain a modified 2 — > 2 part of the total cross section. Note that this convolution also 
includes soft and collinear photonic corrections to the Born/one- loop interference. The complete 
result again contains in addition the 2 — > 3 part, 

(^tot,isn(s, m 2 e ) = o- v+SjIS R+ / dr 3 \M 2 ^ 3 {s)\ 2 . (17) 

J AE-,,A8-i 

3.3 Simulation 

As can be verified in Fig. the resummation approach does eliminate the problem of negative 
weights: shifting the energy cutoff below the experimental resolution, e.g., AE y = 0.5 GeV, such 
that photons are explicitly generated whenever they can be resolved, the subtracted effective 
squared 2 — > 2 matrix element is still positive semidefinite in the whole phase space. Since 
neither the inclusion of the ISR structure function nor the addition of the 2—^3 part introduces 
further sources of negative weights, unweighting of generated events is now possible, so this 
method allows for realistic simulation at NLO. 

After resummation, the only potentially remaining source of negative event weights is the 
soft-noncollinear region. Negative weights are absent as long as 

0(1) x -ln^£ln(A0 7 ) 

7T S 

stays smaller than one, where the 0(1) prefactor depends on the specific process. For the cutoff 
and parameter ranges we are considering here, this condition is fulfilled. 
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A final possible improvement is to also convolute the 2 — 3 part with the ISR structure 
function, 



This introduces another set of higher-order corrections, namely those where after an arbitrary 
number of collinear photons, one hard non-collinear photon is emitted. This additional resum- 
mation does not double-count. It catches logarithmic higher-order contributions where ordering 
in transverse momentum can be applied. Other, logarithmically subleading contributions are 
missed; this is consistent since the genuine second-order part is not calculated anyway. 

4 Results 

4.1 Choosing Cutoffs 

In the kinematical ranges below the soft and collinear cutoffs, several approximations are made. 
In particular, the method neglegts contributions proportional to positive powers of AE^ and 
A# 7 , so the cutoffs must not be increased into the region where these effects could become 
important. On the other hand, decreasing cutoffs too much we can enter a region where the 
limited machine precision induces numerical instabilities. Therefore, we have to check the 
dependence of the total cross section as calculated by adding all pieces and identify parameter 
ranges for AE 1 and A# 7 where the result is stable but does not depend significantly on the 
cutoff values. 

Energy cutoff dependence 

In Fig. |3] we compare the numerical results obtained using the semianalytic fixed-order calcu- 
lation with our Monte-Carlo integration in the fixed-order and in the resummation schemes, 
respectively. Throughout this section, we set the process energy to ^/s = 1 TeV and refer to 
the SUSY parameter point SPSla'. All 2^2 and 2^3 contributions are included, so the 
results would be cutoff-independent if there were no approximations involved. 

The fixed-order Monte-Carlo result agrees with the semianalytic result, as it should be the 
case, as long as the cutoff is greater than a few GeV. For smaller cutoff values the Monte-Carlo 
result drastically departs from the semianalytic one because the virtual correction exceeds the 
LO term there, and therefore the 2^2 effective squared matrix element becomes negative 
in part of phase space. For the Monte-Carlo approach, aiming at unweighting events, the 
integrand is set to zero in regions where it is actually negative, and the result overshoots when 
this happens. 




(18) 
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Figure 5: Total cross section dependence on the energy cutoff AE^ using different calcula- 
tion^ methods: 'sa' (dotted) = fixed- order semianalytic result \15$ using FeynArts/ForrnCalc; 
'fix' (dashed) = fixed-order Monte-Carlo result |2J) using WHIZARD; 'res' (long-dashed) = ISR- 
resummed Monte-Carlo result MM) using WHIZARD; (dash-dotted) = same but resummation ap- 
plied only to the 2—^2 part \1T\) . Statistical Monte-Carlo integration errors are shown. For the 
Monte-Carlo results, the collinear cutoff has been fixed to A0 7 = 1°. The Born cross section is 
indicated by the dotted horizontal line. 



The semianalytic fixed-order result is not exactly cutoff-independent, as one could naively 
expect. Instead, it exhibits a slight rise of the calculated cross section with increasing cutoff; 
for AE^ = 1 GeV (10 GeV) the shift is about 2 %o (5 %o) of the total cross section, respectively. 
While for cutoff values approaching the process energy the soft approximation is expected to 
fail, the rise at small cutoff values is due to the fact that in the soft-photon factor fsoiti—^ 1 ) the 
kinematics is slightly tweaked (necessary to cancel the photon-mass dependence in the virtual 
part), and the cancellation of the logarithmic singularity in AE^ is therefore not exact. The 
mismatch is parameterically of the order AE^/y/s multiplied by the virtual correction; given 
the fact that the virtual correction exceeds the LO term at about AE^ = 1 GeV, we expect an 
error of up to a few %o at that point. For E < 0.1 GeV, this error becomes truly negligible. 

The fully resummed result (118)1 shows an increase of about 5 %o of the total cross section 
with respect to the fixed-order result which stays roughly constant until AE 1 > 10 GeV where 
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Figure 6: Total cross section dependence on the collinear cutoff AO 1 using three different calcu- 
lational methods: sa (dotted) = fixed- order semianalytic result ti,5]/ using FeynArts/ForrnCalc; 
fix (dashed) = fixed-order Monte-Carlo result (0j using WHIZARD; res (long-dashed) = third- 
order ISR-resummed Monte-Carlo result Alfy) using WHIZARD; (dash-dotted) = same, but using 
the first- order ISR-resummed structure function instead. Statistical Monte- Carlo integration 
errors are shown. The soft cutoff has been fixed to AE^ = 10 GeV. 



the soft approximation breaks down. This increase is a real effect; it is due to higher-order 
photon radiation that is absent from the fixed-order calculation. Comparing with the curve for 
2^2 resummation ()17)) . we observe that for AE^ > 1 GeV these higher-order contributions 
are caught by ISR resummation of the 2 — > 2 part, but are transferred to the 2 — > 3 part if the 
cutoff is lowered further, i.e., one radiated photon is resolved. 

At the ILC, with a cross section of more than 100 fb and an integrated luminosity of 1 ab _1 , 
a statistical la fluctuation level of 2.5 %o is reached. Although systematical uncertainties in the 
analysis are likely to be relevant as well, the theoretical prediction of the Monte-Carlo generator 
should aim at matching that precision. To get rid of artefacts of the soft approximation at 
the level of 2 %o, we have to choose AE^ < 0.5 GeV. (For cutoffs lower than 0.1 GeV double- 
precision numerics breaks down, although we still could obtain results by switching to quadruple 
precision.) Resumming photons leads to an increase of 5 %o. 

A common practice is to just convolute the Born part with the ISR structure function, 
leaving all other contributions at fixed 0(a). Our results show that this is insufficient if the 
achievable accuracy is to be exploited. 

To conclude, for the resummation method indicated by Eq. (fTHjl the desired low cutoff 
values are actually acceptable, and the systematic errors induced by the soft and collinear 
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approximations can be suppressed down to the expected level of statistical fluctuations. In 
principle, NNLO corrections and higher order effects of running couplings should be studied for 
a final verdict on the theoretical accuracy. However, we do not expect these corrections to be 
significant; in particular, at ILC energies electroweak Sudakov logarithms are still sufficiently 
small |2Zj • Off-shell and finite-width effects can be taken into account by interfacing the results 
obtained here with the multi-particle event generators presented in |12j . 

Collinear cutoff dependence 

The collinear cutoff A# 7 separates the region where, in the collinear approximation, higher- 
order radiation is resummed from the region where only a single photon is included, but treated 
exactly. We show the dependence of the result on this cutoff is in Fig. 

The plot shows that the main higher-order effect is associated with photon emission angles 
below 0.1°. Cutoff values between 0.1° and 10° are essentially equivalent. To achieve this 
cutoff-independence, collinear terms have to be included in the structure function beyond first 
order (up to third order in our case); using the first-order ISR-resummed structure function 
instead would miss some radiation at low angles # 7 < 1°, cf. the small difference at A# 7 = 0.1° 
between the first- and third-order results in Fig. El For 6* 7 > 10°, the collinear approximation 
breaks down. 

Photon mass dependence 

The dependence on the fictitious photon mass A is eliminated by implementing the soft-photon 
factor / so f t before any further manipulations are done. Therefore, while the photon mass remains 
a parameter in the matrix element code, the result does not numerically depend on it, regardless 
which method has been chosen. 

4.2 Energy Dependence of the Total Cross Section 

Fixing the cutoffs to reasonable values, we can use the integration part of the Monte-Carlo 
generator to evaluate the total cross section at NLO for various energies. This is shown in 
Fig. where we display the LO result together with the NLO result for the fixed-order and 
resummed approach indicated above, respectively. Near the cross-section maximum, the relative 
correction in the fixed-order (resummed) approach is about —5%, approaching —2% (—1.5%) 
at y/s = 1 TeV, respectively. Near threshold and at asymptotic energies, the relative NLO 
correction is larger in magnitude. 

4.3 Simulated Distributions 

The strength of the Monte-Carlo method lies not in the ability to calculate total cross sec- 
tions, but to simulate physical event samples. We have used the WHIZARD event generator 
augmented by the effective matrix elements and structure functions as introduced above, to 
generate unweighted event samples for chargino production. 
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Figure 7: Left: o~ tot in fb as a function of \fs; right: relative correction to the Born term 
in percent. LO (dotted) = Born cross section without ISR; NLO/fix (dashed) = fixed-order 
approach; NLO/res (full) = resummation approach. We choose AE^ = 0.005-y/i and A# 7 = 1° 
as cutoffs separating the hard (non-collinear) from the soft (collinear) regions, respectively. 
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Figure 8: Polar scattering angle distribution for an integrated luminosity of 1 ab~ at yfs = 
1 TeV. Left: total number of events per bin; right: difference w.r.t. the Born distribution. LO 
(dotted) = Born cross section without ISR; fix (dashed) = fixed-order approach; res (full) = 
resummation approach. Cutoffs: AE^ = 3 GeV and A6* 7 = 1°. 
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To evaluate the importance of the NLO improvement, in Fig. |H1 we show the binned dis- 
tribution of the chargino production angle as obtained from a sample of unweighted events 
corresponding to 1 ab _1 of integrated luminosity. The collider energy has been set to 1 TeV, 
the SUSY parameter point SPSla' is the same one as for the previous plots. With cutoffs 
A# 7 = 1° and AEj = 3 GeV we are not far from the expected experimental resolution, while 
for the fixed-order approach negative event weights do not yet pose a problem. (As discussed 
above, for the resummation approach decreasing cutoffs further is possible and preferred, but 
choosing lower values would invalidate the fixed-order approach for the comparison.) 

The histograms illustrate the fact that NLO corrections in chargino production are not 
just detectable, but rather important for an accurate prediction, given the high ILC luminosity. 
The correction cannot be approximated by a constant i^-factor but takes a different shape than 
the LO distribution. The correction is positive in the forward and backward directions, but 
negative in the central region. The effects of photon resummation are not as striking, but still 
statistically significant; they are visible mostly in the central-to-forward region. Apparently, to 
carefully choose the resummation method and cutoffs will be critical for a truly precise analysis 
of real ILC data. 

5 Summary and Outlook 

We have presented results obtained from implementing NLO corrections into a Monte-Carlo 
event generator for chargino pair-production at the ILC. On top of the genuine SUSY/electroweak 
corrections, we have considered several approaches of including photon radiation, where a strict 
fixed-order approach allows for comparison and consistency checks with published semianalytic 
results in the literature, while a version with soft- and hard-collinear resummation of photon 
radiation not just improves the numerical result, but actually is more straightforward to im- 
plement and does not suffer from negative event weights in or near the accessible part of phase 
space. 

A careful analysis of the dependence on the technical cutoffs on photon energy and angle 
that slice phase space in regions, reveals uncertainties related to higher-order radiation and 
breakdown of the soft or collinear approximations. For the level of precision required by ILC 
analyses, the cutoffs have to be chosen rather low. Resummation of photon radiation is required 
not just for precision, but also to get rid of negative event weights in the simulation. In the 
Monte-Carlo event generator WHIZARD we have thus implemented the NLO result with higher- 
order resummation in critical regions. The generator accounts for all yet known higher-order 
effects, allows for cutoffs small enough that soft- and collinear-approximation artefacts are 
negligible, and explicitly generates photons where they can be resolved experimentally. 

The generator that we have constructed should be regarded as a step towards complete NLO 
simulation of SUSY processes at the ILC. If charginos happen to be metastable, it already pro- 
vides all necessary ingredients. Beam effects (beamstrahlung and energy spread, polarization) 
are available for simulation and can easily be included. However, charginos are metastable only 
for peculiar SUSY parameter points; in general we have to take into account chargino decay 
and the corresponding additional NLO corrections. These we have to match with off-shell and 
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background effects, already available for simulation in WHIZARD. Furthermore, in the threshold 
region the Coulomb singularity calls for resummation, not yet accounted for in the program. 
These lines of improvement will be pursued in future work. 
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A Soft and collinear photon factors 



The soft-photon factor that allows us to eliminate the ficticious photon mass A in favor of a 
physical photon-energy cutoff AE^ is given by the integral 28 22 29J 

d?k (±)p i p j Q i Q j 



soft 



a 
2^ 



|k|<AE, 



2u k (pik)(pjk) 



(19) 



Pi (k) denote the electron/chargino (photon) four- vectors, respectively, while 

uj k = Vk 2 + A 2 



(20) 

is the energy of a photon regularized by the photon mass A. Qi are the corresponding charges. 
The sign is — for incoming and + for outgoing particles, respectively. 

For conserved helicity, the collinear radiation of one photon can be approximated by con- 
voluting the no-photon matrix element with the structure function (see, e.g., [80] ) 



t+f \ i] l+x 2 
f W = T 



4 1 



x 



(21) 



Here, with h^ ax = p° e A# 7 being the collinear cutoff, the expansion parameter 77 is defined as 



2a 

v = — 



7T 



hi 



Am] 



{AO, 



The helicity-flip structure function, 



/-(*) 



a 



[22) 



(23) 



does not contain a logarithmic enhancement with A# 7 , so the helicity-flip part of collinear 
radiation is subdominant. 

In the soft-collinear region x « 1, 7 ~ 0, the leading logarithms in A9 1 (i.e., powers of 77) 
can be resummed to all orders |2H]. Matching this with the complete x-dependent expressions 
for hard collinear radiation to first, second, and third order in rj, Skrzypek and Jadach obtained 
the ISR structure function 26 
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(24) 
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B Two-Photon Phase-Space Regimes in the Resumma- 
tion Method 



The class of processes we are considering exhibits the usual infrared and collinear (for m e = 0) 
singularities in individual contributions to the physical result, which we treat by standard 
methods. In the fixed-order approach where only one photon is present, the consequences of 
this phase-space slicing are evident: (i) In the soft region, the photon energy is set to zero 
in the matrix element, and the real contribution is cancelled against the IR-divergent part of 
the virtual correction, neglecting corrections proportional to the photon energy, (ii) In the 
hard-collinear region, the photon-emission part of the matrix element is replaced by a structure 
function, neglecting corrections proportional to the separation angle. As long as the problem 
of negative event weights can be ignored, and numerics is not an issue, the description at fixed 
NLO always improves if any of those cutoffs are lowered. 

However, the resummation method described in Sec. El involves all orders of soft-photon 
radiation. Shifting cutoffs changes the type of higher-order contributions that are included, so 
lowering cutoffs not necessarily improves the description. 

To clarify this issue, let us focus on the 0(a 2 ) correction, i.e., two photons (real or virtual). 
Since we do not consider two-loop diagrams, this correction is not completely accounted for, 
but dominant contributions are included. 

In the resummation method, there are three different ways of dealing with real and virtual 
photons: 

(a) soft approximation j2I|: 

describes collinear and non-collinear soft photons; neglects contributions cx is com- 

bined in the sequel with the soft photonic part of the one-loop matrix element 

(b) ISR |2E|: 

describes collinear real and virtual photons; neglects interference terms in photon emis- 
sions. Assumes fc T -ordering of the emitted photons, i.e. for j > i: kj > kj , and in n th 
order: YJl=x kj < k^ ax , where k^ ax is fixed. 

(c) real emission given by exact (hard non-collinear) matrix element A^2^3 

Considering now the treatment of two photons (i.e. 0(a 2 ) corrections), at least one of 
the photons is always described by the ISR structure function. But when the Born term is 
convoluted with the ISR function, there are also two-photon contributions described solely by 
the ISR. We have to distinguish between the cases where (i) the two photons are attached to 
the same or (ii) to different incoming particles. In case (i), we consider the three terms 

0(a 2 ) lSR - 0{a) ISK 0(a)^ + 0(a) 1SR 0(a) so{t . (25) 

The first term contains all pairs of collinear photons from the ISR, /c T -ordered; the last term 
contains a first photon from ISR and a second one from the soft-photon factor (SPF, which in 
the following is understood to include the soft photonic one- loop contribution). The term in the 
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middle is the subtraction to avoid double-counting of soft photons described in subsection 13.21 
Here both photons are from the ISR, the first one with arbitrary energy, the second one soft. 

If the second of the considered photons is soft, and both are fc T -ordered, then there is an 
exact cancellation between the first two terms. For non £; T -ordered photons, the first term gives 
no contribution, and there is a cancellation between the second and third term, which results 
in a difference between the "exact" SPF expression and the ISR LLA term for the incoming 
particle 



Since here there are always two different structure functions involved, k -ordering is absent, 
and after a cancellation of soft terms one is left with 



which is up to the missing terms AiA 2 equivalent to an SPF description for both legs. 

We now investigate the changes induced by raising one of the two cutoffs. Generally, for the 
contributions with two real photons, we loose contributions if we lower the cutoffs since double 
photon radiation is not accounted for by the 2 — > 3 matrix element. Convoluting also the 2-^3 
matrix element with the ISR structure function as proposed in Equ. (fTHjl gives contributions 
with a collinear and hard non-collinear photon and cures the problem. But still, raising the 
cutoff A6* 7 opens up phase space for the first photon. Thus, raising the cutoffs gives a better 
description of these contributions as long as the collinear approximation is valid. The same 
holds for the energy cutoff, but we see from Figs. El and El that the soft approximation fails 
much earlier than the collinear description. 

By including the next order of real photon radiation explicitly, i.e., 



we can further improve the description of this part of two-photon phase space. This contri- 
bution, which is however tiny for the cutoff values considered here, can easily be added using 
WHIZARD tree-level event generator. 

For completeness, we finally discuss the reshuffling of contributions in the overlap region 
of the soft-collinear and hard-collinear (ISR) descriptions. If we raise AE^ while keeping A# 7 
fixed, photons that have been hard now become soft. In the case (ii) (photons radiated from 
two different external particles), a photon which has been described by the structure function, 
comes now with the SPF. For the case (i), we have to distinguish whether the two photons 
are fc T -ordered or not. If they are, the description again changes from the ISR to the SPF. If 



Aj = 0(a)j, Boft - O^j^sR. 



In the case (ii), we write the terms schematically as 



C(a)i,iSRC ) (a)2,isR 

+ ©(aKisR (C(a) 2)S oft - 0(a) s 2 % R ) 



(26) 



AiO(a) 2 ,i S R + C(a)i,i S RA 2 + 0{a)x^O (a) 2jISR , 
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there is no k -ordering, then the photons either change from hard+soft to soft+soft, which is a 
smooth transition where only the last two terms of Equ. (J23j) are involved, or the second photon 
changes to soft for the combinations hard+hard or soft+hard. In that case there appear new 
contributions of the form AO(a)isR, which have not been there before. 

Raising A6 1 , while keeping AE 1 fixed, shuffles photons from a non-collinear to a collinear 
description. The interesting region is for photons that lie in the soft regime near the limit of the 
soft-collinear regime and change into the latter after raising the angle cutoff. For A; T -ordered 
photons the cross-over is smooth and the second photon stays with the SPF, while for non- 
ordered photons the description switches from the SPF to the the difference between ISR and 
SPF. For the case (ii) of radiation from different legs, the description always remains with the 
SPF. 

For more details see j31j . 



C SPSla' 

The SUSY parameter point SPSla' is defined in Ref. [3J; it is a SUGRA-type scenario derived 
from the parameter set 

m = 70 GeV, m 1/2 = 250 GeV, tan (3 = 10, 

\i > 0, A = -300 GeV. (27) 

The precise spectrum and coupling parameters are computed using the renormalization-group 
evolution code of Ref. the values can be found in Ref. [3]. The chargino masses and widths 
are also listed in Tabled 
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